Vector chiral order in frustrated spin chains 
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By means of a numerical analysis using a non-Abelian symmetry realization of the density matrix 
renormalization group, we study the behavior of vector chirality correlations in isotropic frustrated 
chains of spin 5 = 1 and S = 1/2, subject to a strong external magnetic field. It is shown that 
the field induces a phase with spontaneously broken chiral symmetry, in line with earlier theoretical 
predictions. We present results on the field dependence of the order parameter and the critical 
exponents. 
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I. INTRODUCTION 

The so-called vector chirality in quantum spin chains is 
defined as the vector product of two adjacent spins along 
the chain: 

In a chirally ordered state spins have a tendency to "ro- 
tate" in some preferred plane in a certain preferred di- 
rection (clockwise or counterclockwise). Chiral phases 
in quantum spin chains were predicted long ago!* 2 - and 
have attracted considerable interest recently 3 -^^! 7 ^ af- 
ter they have been found numerically in frustrated chains 
with easy-plane anisotropy.^ 7 - As noted by Villain,— the 
chiral order should survive at finite temperature in the 
presence of three-dimensional interactions without tran- 
forming into a usual helical long-range order: at finite 
temperatures the chirality correlation length is much 
larger than the spin correlation length, so with decreasing 
temperature chiral order should set in before spin order 
does; there are experimental indications that the chiral 
order may exist in the quasi-one-dimensional anisotropic 
organic magnet Gd(hfac)3NITiPr. 9 The projection of the 
vector chirality k on the direction of the applied field 
could be experimentally detected with the help of polar- 
ized neutrons.— 

In all known cases of numerically confirmed existence 
of chiral states, the preferred plane for spin rotation is 
chosen by some anisotropy of the easy-plane type, and 
the chiral phase disappears in the isotropic limit.— Re- 
cently, it has been predicted 8 that in isotropic frustrated 
chains the chiral phase may appear in presence of an 
external magnetic field, strong enough to close the spec- 
tral gap. In such a state, the system approximately de- 
couples into a gapped antisymmetric sector and a gap- 
less symmetric sector, the latter being described by the 
Tomonaga-Luttinger liquid (TLL). An alternative two- 
component TLL scenario 11 ' 12 ! 13 assumes the existence of 
the Tomonaga-Luttinger liquid in both sectors and im- 
plies absence of the chiral order. 

The phase diagrams of the antiferromagnetic zigzag 
spin chains in applied field have been studied numeri- 
cally in Ref. [U for 5 = 1/2, and in Ref. for 5 = 1. 



However, both works focused only on the magnetization 
process and did not check the presence of the chiral or- 
der. The theoretical analysis of Ref. H involves some un- 
controlled approximations (mean-field decoupling of the 
"twist" term) and, to our knowledge, the chiral order 
in isotropic spin chains has never been directly probed 
numerically (the only exception being the calculation of 
short-range correlations in a 5 — 1/2 chain^), so the 
question of the correct scenario remains unsettled. 

Several materials are known which realize isotropic 
zigzag spin chains^ among them, (N2H5)CuCi3 can be 
viewed as a promising candidate for experimental stud- 
ies, since its small exchange constants make it feasible to 
reach magnetic fields comparable to the gap. 

In this paper we present a study of vector chirality 
correlations in the isotropic 5 = 1 and 5=1/2 zigzag 
chains in the presence of applied magnetic field, using a 
powerful non-Abelian symmetry realization 17 of the den- 
sity matrix renormalization group technique 18,19 in its 
matrix product state formulation.— It is demonstrated 
that the chiral order does exist in the high-field phase, 
both for 5=1 and 5 = 1/2, and the behavior of chiral 
correlations is in a qualitative agreement with the ex- 
pectations following from the theoretical analysis of Ref. 
8. This implies that a chiral one-component Tomonaga- 
Luttinger liquid scenario is realized. 



II. THEORETICAL ESTIMATES 

We consider the model of a frustrated antiferromag- 
netic spin chain, defined by the Hamiltonian: 

H = J\ ^ S n ■ S n+ i + J2 ^ S„ ■ S n+ 2 — H ^ S„ (1) 



where S n are spin-5 operators at the n-th site, J\ > 
and J2 > are the nearest and next-nearest neighbor 
exchange constants, respectively, and H is the external 
magnetic field, assumed to be applied along the z axis. 

In case of 5 = 1, at H = the ground state is al- 
ways gapped: for small frustration parameter a = J2/J1 
one remains in the Haldane phase characterized by the 
long-range string order, while for a > a c « 0.75 there 
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FIG. 1: A zigzag spin chain and the notation adopted in 



is a first-order transition into another gapped state, the 
so-called "double-Haldane" (DH) phase where the string 
order disappears with a finite jump, giving way to a more 
complicated hidden order i2I When the applied field ex- 
ceeds the critical value H = H c (H c is obviously equal to 
the gap at H = 0), the system acquires finite magnetiza- 
tion. There is another special field value, the saturation 
field H s above which the spins are fully polarized. In the 
5 = 1/2 case, the ground state at H — is gapless for 
a < a c ~ 0.24, and a > a c corresponds to a gapped 
dimerized phase. 

We are interested in the properties of the partially 
magnetized state in the field range H c < H < H s , and we 
assume that a is large enough to make the system gapped 
at H = 0. In the limit a>l the system may be viewed 
as two weakly coupled chains. A single S — 1/2 chain in 
external field has been extensively studied ) 22 i 23 ' 24 as well 
as its S = 1 counterpart i 25 i 26 ' 27 i 28 Above the first critical 
field H c i the low-energy physics of a single chain is well 
described in terms of the effective Tomonaga-Luttinger 
liquid (TLL) theory, described by the Hamiltonian 

Htl\BM = \ J dx{^{d x <j>) 2 + K(d x 6) 2 }. (2) 

Here K is the so-called TLL parameter, v oc J2 is the 
Fermi velocity, (f> is the bosonic field (compactified by 
the condition (p = if + and 9 is its dual satisfying 
the commutation relations [4>(x), 9(y)\ — iQ(y—x), where 
0(x) is the Heaviside function. 

In the continuum limit, the lattice spin operators both 
for S = 1/2 and S — 1 can be represented^ 2 ^ through 
the bosonic field ip and its dual 9: 

S z a (x a ) = M + -?=d x <j> a (x a ) 

+ A 3 sin {irMx a + V^M x a) } + (■■■) (3) 
S+(x a ) = e i™/2 e iV*e«(* a ){ Al 

+ A 2 sin {-KMx a + V^K(j) a {x a )) } + (■■■), 

Here a = 1,2 labels the two chains, the space coordinate 
x is defined along the zigzag path as shown in Fig.[T] the 
lattice sites correspond to X\ t 2 = x =F i, M is the ground 
state magnetization per spin in units of saturation (which 
corresponds to the filling factor in the TLL model), and 
Ai are nonuniversal amplitudes. In the case of S = 1 
the dots denote additional operators which correspond to 
massive fields connected to the high-energy S z — 0,-1 
magnon branches of the Haldane chains. 



The theory parameters M, K, and v should be un- 
derstood as functions of the field H; for S = 1 their 
behavior can be extracted from the comparison of the 
TLL theory predictions with the numerical resultS ) 26 i 29 
while for S = 1/2 it is available from the exact Bethe 
ansatz slution^ 2 - 

For S = 1 the most important featur o 25 i 26 ' 29 is that 
generally the TLL parameter K > 1 at H > H c \, 
non-monotonically depends on H and tends to the free 
fermion value K = 1 both at the first critical field 
H = H c i and at the saturation field H = H s i. In con- 
trast to that, for S — 1/2 the TLL parameter K < 1 
and is a monotonically increasing function of the applied 
field^^ 

For the description of weakly coupled chains it is con- 
venient to introduce the symmetric and antisymmetric 
combinations of the bosonic fields 



= (0i ± ' 



The longitudinal (S Z S Z ) part of the zigzag exchange, 
apart from producing terms of the type (d x fi)(d x f2) 
which lead to a splitting of the TLL parameter values for 
the symmetric and antisymmetric sectors, 



K±ttK[l±2K/(nva))] 1/2 , 
v± « v[l±2K/(Trva))] 1/2 , 



(4) 



yields another contribution proportional to cos [^/8Trip-~ 
7rM] . The scaling dimension of this latter perturbation 
is 2K_ . In the case of S = 1 it is irrelevant since K > 1 
and so can be neglected; in contrast to that, for S = h 
chain K < 1, and this operator is a relevant perturbation. 
Thus, as pointed out in Ref. [^, for S = 1 the leading con- 
tribution to the interaction is given by the "twist term" 
produced by the transversal part of the zigzag exchange. 
For S =1/2 the twist term competes with a relevant op- 
erator and can only win if K_ (H) is above a certain criti- 
cal value K c ; for K < K c the so-called "even-odd" phase 
is realized, whose dominant correlations are of the spin- 
nematic (or XY2 in the nomenclature of Ref. 30) type. 
A characteristic feature of the "even-odd" phased are 
the AS Z = 2 steps in the magnetization curve S* ot (H). 
The "even-odd" phase has been also observed^ in zigzag 
chains with ferromagnetic nearest-neighbor exchange. 

In the chirally ordered phase the twist term is the most 
relevant perturbation, so one obtains the same effective 
Hamiltonian for S = 1 as well as for S = 1/2: 



HcS = ^ Ti-TL [(fig , 9 g] + Hj, 



<j=± 



TLw 



dxsin (v2ir9-)(d x 9- l 



(5) 



Mean-field decoupling of the twist term in the spirit of 
Ref. y then leads to the conclusion— that both (d x 9 + ) and 
(sin ( k \/2tt9J)) become nonzero, and the antisymmetric 
sector becomes gapped. 
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FIG. 2: Numerically calculated magnetization curve S* ot (H) 
for a S = 1 zigzag chain of length L = 128 with frustration 
parameter a = J2/J1 = 1. 



One should mention that the above description makes 
sense only when we are far enough from the critical fields 
H c or H s : the theory is applicable only up to the en- 
ergies of the order of the bandwidth v, and v — > if 
H — > H C ,H S . The formulae Q indicate that the system 
becomes unstable against phase separation^ as soon as 
v < 2K/(na). 

The components of the chirality operator n can be ex- 
pressed through bosonic fields. The longitudinal part of 
the chirality can be obtained in the following form: 



k z (x) = sin(\/27r0_){A? 



A 2 -i 
+ -± cos (V&Kip + + 27rMx) I + (•••), (6) 

where dots denote massive fields (the most impor- 
tant contribution of that sort is proportional to 
(— 1)* cos Cs/2n0— ) {d x 6+) ) and operators with higher 
scaling dimensions. The leading contribution to the long- 
distance correlator is thus given by 



MO 2 
4 

2nMx) 



C 2 cos(27rMa;) 



(7) 

where £ is the largest correlation length determined by 
the gap in the antisymmetric sector. (For S = 1 there 
is also another, much smaller, characteristic correlation 
length £ which is determined by the high-lying excitation 
branches that are neglected in the bosonization formulae 
©; it roughly corresponds to the correlation length of the 
Haldane chain at zero field, typically a few lattice con- 
stants.) For 5=1, although K > 1 for H c < H < H s , 
the parameter K + , according to (U]), is renormalized to 
smaller values when the zigzag coupling is switched on, so 
the two decaying contributions in {7} may be competing 
with each other. In S = 1/2 chain the oscillating contri- 
bution always has the slowest decay since K + < K < 1. 
In a similar way one can obtain the transversal chirality 
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FIG. 3: (Color online) Typical DMRG results for the chirality 
correlator (|13[) and results of fitting it to the form (|14[) . The 
error bars 3 — (here smaller than the symbol size) indicate the 
variance of F(x — n — n') calculated from averaging over the 
initial and final points n, n' . 



component 

K+(x) = 2AlMsm(^/^/2e-) 

x exp[— +l ^-0+} + (•••)• (8) 

It is easy to see that the leading term in n + (x) is simply 
proportional to — St- The leading contribution to the 
corresponding asymptotic correlator is slowly decaying 

M 2 

(k+(x)k-(0)) ex Aj xl/(4K+) exp{iQx}, x > £, (9) 
and incommensurate, with the wave vector given by 

The above expressions ([7} and © are expected to be 
valid in the limit a 3> 1, and for H not very close to 
the critical fields H c , H s (in the vicinity of the critical 
field the bosonization approach becomes hardly applica- 
ble since the effective bandwidth goes to zero). Close to 
the saturation field H s , a large- S analysis^ allows map- 
ping the system to an effective model of two bosonic 
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species with repulsive interaction, which condense driven 
by the magnetic field playing the role of the chemical 
potential. The repulsion turns out to be strong enough 
to satisfy the phase separation condition, so only one of 
the species condenses, the other condensate is depleted, 
so one deals in fact with the one-component pseudo- 
condensate whose physics is again described by a (one- 
component) TLL. The asymptotic form of the longitudi- 
nal chirality correlator for H close to H s has been pre- 
sented in Ref . d: 

(^(zK(0))-^-^, (10) 

with Kg oc (H s — H). The leading contribution to the 
transversal chirality is proportional to the bosonic field it- 
self, so its correlator takes the following asymptotic form: 

< K +(*) K -(0))^^g^e^. (11) 

Here K' is another TLL parameter, the characteristic 
wave vector Q' is given by the expression for the pitch of 
the classical helical state: 

Q' = ±(tt — arccos(l/4a)), (12) 

and the amplitude C oc (i^-iT) 1 / 2 " 1 /^'). As the field 
approaches the saturation point, H — » H s , the value of 
K' tends to 1, so the amplitude C vanishes. 

III. RESULTS OF NUMERICAL ANALYSIS 

We have studied the S = 1 and 5=1/2 zigzag chain 
model given by Jl} using the DMRG method in its matrix 
product state formulation, making full use of the non- 
Abelian SU(2) symmetry. 

For a full description of the DMRG technique, 18 , we re- 
fer the reader to the review Ref. LLiJ As discussed in Ref. 
I20L the formulation in terms of matrix product states is 
very convenient but for the calculation of ground states 
does not lead to substantially better results. The decisive 
poin t 17 i 20 is the use of the non-Abelian symmetry SU (2) 
instead of the Abelian U(l). While the magnetic field 
H breaks SU(2) symmetry, the fact that the Zeeman en- 
ergy term commutes with the rest of the Hamiltonian 
makes it possible to take take the influence of the mag- 
netic field into account by calculating the ground state 
of the model in a sector with the given total spin Stot- 
The advantage of the method lies in a drastic reduction 
of the number of states m which is necessary to describe 
the system, because non-Abelian symmetry allows one to 
calculate using representatives of groups of states of the 
same total spin: essentially, one treats the multiplet of 
states of the same total spin as a single representative 
state. Comparing to the Abelian version of the method 
which only uses the £7(1) symmetry, the improvement in 
efficiency can be several orders of magnitude, depending 
on the problem. For the zero-field groundstate of the 




FIG. 4: (Color online) (a) the square chirality order param- 
eter Kq as a function of magnetization M = Stot/L for a 
S = 1 zigzag chain with a — 1, extracted from fits of chiral- 
ity correlation functions; (b) zoom-in of the same picture in 
the vicinity of M = 1, where finite size effects become impor- 
tant. The inset shows the L scaling of the point M c where 
the long-range chirality order disappears in a finite system of 
size L. 

zigzag chain, the effective improvement in the number of 
states is a factor ~ 3 (for S = 1/2) or ~ 4 (for S = 1), 
leading to a reduction in the computational effort by a 
factor ~ 27 and ~ 64 respectively. The relative efficiency 
decreases as the magnetic field is increased, but even for 
rather high fields the improvement is appreciable. 

A slight disadvantage is that the non-Abelian method 
allows to compute only reduced matrix elements (in the 
sense of the Wigner-Eckart theorem). In our case, since 
the chirality is a vector, the correlator that is by far the 
easiest to calculate is the rotationally invariant scalar 
product, 

F(n-n') = {K(n)-K(n')), (13) 

which is a mixture of the longitudinal and transversal 
contributions^ This obviously makes the analysis of the 
numerical data more difficult: in our case, from the the- 
oretical analysis it follows that the longitudinal chirality 
correlations (k z (x)k z (0)) decay to their asymptotic value 
much faster than the transversal ones. Thus it turns out 
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FIG. 5: (Color online) Behavior of the transversal chirality 
correlations exponent ij for a S = 1 chain with a = 1 as a 
function of magnetization M = Stot/L, extracted from fits of 
the correlation function 1)131) to the functional form (| 14[) . The 
error bars shown correspond to the uncertainties of the fit. 



FIG. 6: (Color online) Behavior of the oscillations amplitude 
A as a function of the magnetization M — Stot / L for a S — 1 
chain with a = 1. One can see that the scaling A oc M 2 
suggested by © is only applicable for small M. 



to be practically impossible to extract the characteristic 
decay exponent ry 2 for the longitudinal chirality correla- 
tions from the F(x) data, and one can only try to esti- 
mate the exponent rj — rj xy of the transversal chirality 
correlations. 



A. S = 1 zigzag chain 

We have studied spin-1 zigzag chains with the frustra- 
tion parameter a = 1, for several chain lengths L rang- 
ing from 64 to 192. For our calculation, even within the 
SU(2) method we needed a relatively large number of 
representative states(from m ~ 400 to to ~ 2000, de- 
pending on L and Stot) to reach good convergence. In 
the U(l) formulation, this corresponds to an m of up to 
8000, making an Abelian calculation much more difficult. 

As one can see from the numerically calculated mag- 
netization curve shown in Fig. [5J at a = 1 the S Z (H) 
dependence is featureless and shows neither plateaux, 
nor cusps, nor AS* 2 = 2 steps characteristic for even-odd 
phase, in accordance with the results of Rcf. 

We have computed the chirality correlator (fl3|) in the 
ground states of a large number of sectors with certain 
total spin quantum number Stot- When computing F{n— 
n'), it was averaged over the starting and final positions 
n, n', and care was taken to stay in the bulk, away from 
the ends of the chain. The DMRG data for the correlator 
has been fitted to the power-law form 

Acos[ q ( X + S)} 

suggested by ([7]) , © . The introduction of a finite phase 
shift 6 is necessary due to the open boundary conditions. 
Typical fits are presented in Fig. [3] From those fits we 



have extracted the behavior of the chirality order pa- 
rameter Kq and the exponent rj as functions of the chain 
magnetization M = Stot/L, shown respectively in Fig. U] 
and Fig. [5] The fitted wave vector q only weakly depends 
on the magnetization: as M changes from to 1, q varies 
from 1.79 to 1.83, which favorably compares to the clas- 
sical value Q' « 1.82 obtained from (fT2|) at a = 1. The 
phase shift 8(a = 1) w 1 ± 0.05 is also practically inde- 
pendent of M . The behavior of the oscillation amplitude 
is shown in Fig. [SJ the scaling A cx M 2 suggested by © 
is indeed observed for small M, and strong deviations 
appear for M > 0.3. 

One can see that for the bulk of M values the order 
parameter Kq has practically reached convergence already 
at L = 128, so there is no need to perform the finite-size 
scaling. The only region where finite-size effects remain 
strong is M — > 1 (see Fig.^J}): the Kq(M) dependence at 
finite L shows kq vanishing at some M = M c ^ 1. The 
finite-size scaling of M c (see the inset of Fig. [4]b) shows 
that there is no trend to convergence even for L = 192, 
although there is a sizeable increase in M c towards 1 with 
increasing the size L. From the theoretical analysis one 
expects M c — > 1 for L — > oo, however, studying this limit 
numerically can be quite difficult: since for M — > 1 one is 
very close to the fully polarized state, the actual number 
of particles in the problem is the number of magnons 
L(l — M), which has to be large enough to observe a 
phase transition. 

The quality of fits deteriorates for very large Stot i since 
the overall scale of the chirality correlations goes to zero 
as M — * 1. The fits become less reliable for small Stot as 
well, for the following reason: with (jT5J) one attempts to 
fit just the oscillating (transversal) part of the chirality 
correlator. For small Stot the amplitude of the oscillating 
part is small, while the gap in the antisymmetric sector 
is small, so the oscillations appear on top of the exponen- 
tial decay characterized by two very different correlation 
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FIG. 7: Typical DMRG results for the cliirality correlator 
(|13l) of a S — 1/2 zigzag chain with a = 1. The error bars 
have the same meaning as in Fig. [3] (a) (Color online) A 
point in the chiral phase. The solid line shows a fit to (I14|l ; 
(b) a correlator in the nonchiral "even-odd" phase. 
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FIG. 8: (Color online) Behavior of the square chirality k§ 
as the function of magnetization M = 2Stot/L for a S = 
1/2 zigzag chain with a = 1, extracted from fits of chiral 
correlation functions. 
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lengths £ (for a = 1 one has £ « 7) and £ ^> £. It thus 
becomes a numerically ill-posed problem to filter out the 
power-law decaying oscillating part on top of such a back- 
ground. 

The extracted critical exponent slowly changes with M 
and lies in the range rj ~ 0.3 -j- 0.6, which qualitatively 
agrees with the theoretical estimates predicting that for 
a ^> 1 it should vary from approximately 0.25 to 0.5; the 
error bars shown in Fig. [5]) are in fact only of indicative 
nature since they only show uncertainties of the fit to 
the fixed fit function and do not take into account the 
variations of the fit parameters which would result from 
adding subleading (e.g., exponentially decaying) contri- 
butions to Eq. (Til)) . 

The chosen value of a = 1 is rather small and does 
not allow a direct comparison of r] with the theoreti- 
cal value (4K+) -1 : if one naively tries to extract the 
velocity parameter v assuming that K + is given by ^ 
and using the data of Ref. [2!^ for the K(M) dependence, 
the obtained values of v fall below the phase separation 
threshold v c = 2K/ (ira) . We have refrained from study- 
ing chains with a>l since, on the one hand, increasing 
a causes a dramatic increase in numerical effort, and on 
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FIG. 9: (Color online) Behavior of the transversal chirality 
correlations exponent 77 for a S = 1/2 chain with a = 1 as the 
function of magnetization M — 2Stot/L, extracted from fits 
of the correlation function (1131) to the functional form (|14|) . 
The error bars shown correspond to the uncertainties of the 
fit. 



the other hand, the chirality correlators become more and 
more "polluted" by the exponentially decaying contribu- 
tions with ever larger correlation length £. 



B. 5 = 1/2 zigzag chain 

We have also computed the chirality correlation func- 
tion (TrU|) for 5=1/2 zigzag chains of length L = 168, 
256, with frustration parameter a = 1. Typically, 300 
to 400 representative SU (2) states were kept in the cal- 
culation. The magnetization curve for a S = 1/2 zigzag 
chain with a = 1 has been presented in Fig. 3a of Ref. 
[ill According to the phase diagram obtained in Ref. 
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Il3l . at a — 1 the 5 = 1/2 chain exhibits several phases 
with varying the applied field H. Some of those phases, 
namely, the "even-odd" phase and the plateau phase, ob- 
viously do not possess any chiral order. "Suspect" with 
respect to chirality are only two regions marked "TL2" 
in Ref . HH and identified as a two-component Tomonaga- 
Luttinger liquid phase. Indeed, we observe a finite value 
of vector chiral order in both TL2 regions (although not 
in the entire high-field TL2 piece, see below). 

A typical example of the correlator in the chiral phase 
is shown in Fig. [7^; for a comparison, we also show a 
correlator in the "even-odd" phase which is nonchiral. 
The maximal magnitude of chirality in 5 = 1/2 chain 
is roughly one order of magnitude smaller than in the 
5 = 1 case. To analyse the chirality correlators, we had 
to use relatively large chain lengths, because the gaps are 
smaller than in the S = 1 case, and the results for small 
L are polluted by slow exponentially decaying contribu- 
tions. 

We have employed the same fitting procedure as de- 
scibcd above for the 5 = 1 chain, and analyzed the 
behavior of the chirality order parameter k,q and the 
critical exponent rj as functions of the magnetization 
M = 2S to t/L. The results are shown respectively in Fig. 
|8]and Fig. [9] For the low-M chiral region, the amplitude 
of oscillations in the correlation function turns out to be 
too small to extract the exponent 77 with any reasonable 
accuracy, so for that region we were only able to extract 
the order parameter Kq. 

The boundaries of the low-field piece of the chiral phase 
coincide with the low- field TL2 region of Ref.[l3|. Surpris- 
ingly, this is not the case for the high-field piece: while 
its lower boundary reasonably agrees with the transition 
from "even-odd" phase to TL2, its upper boundary lies at 
a finite M = M c ~ 0.75 and not at M — 1 as one expects 
from the theoretical analysis. It is worth mentioning that 
the magnetization curve of the 5=1/2 chain at a = 1 
(see Fig. 3a of Ref. [13|) seems to exhibit a weak feature 



around M ~ 0.75, namely a fast growth of the second 
derivative d 2 M/dH 2 . At M — > M c the critical exponent 
r) tends to 1/2, the value which is expected theoretically 
close to the saturation field. 

In contrast to 5 = 1, where the respective boundary 
exhibited strong finite-size scaling, in the 5 = 1/2 case 
we have not observed any significant change of M c with 
increasing L from 168 to 256, as seen from Fig. [8] We 
have found no chiral order for L = 516 chain with 5t t = 
205, which means that even for such a long chain M C (L — 
516) < 0.787. On the basis of available data, one can 
conclude that the 5=1/2 chain might possess another 
nonchiral phase close to the saturation field. The nature 
of this phase needs further investigation. 



IV. SUMMARY 

We have studied spin-1 and spin-1/2 isotropic anti- 
ferromagnetic zigzag chains in strong magnetic fields by 
means of the matrix product density matrix renormal- 
ization group technique. Existence of a phase with field- 
induced vector chiral order is established for 5 = 1 as well 
as for 5 = 1/2, and the behavior of the order parameter 
and its correlations as functions of the magnetization is 
analyzed. The chiral phase is gapless and corresponds 
to a one-component Luttinger liquid, thereby confirming 
the scenario proposed in Ref. [8|. 
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In that way, the error bars shown do not reflect the actual 
error of the DMRG calculation, which is much smaller, but 
merely reflect the strength of finite-size effects. 
By splitting the vector operator into components that 
change the total spin by AS = 0, ±1, the components of 
the correlator transverse and longitudinal to the field can 
be measured using only rotationally invariant scalar op- 
erators, which are expressible in terms of reduced matrix 
elements of irreducible tensors, but in this form they are 
quite complicated and difficult to construct. 



